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Abstract 

Tam and Auriault 5 successfully predicted the acoustic spectrum at 90° to the axis of a 
high speed air jet by using an acoustic equation derived from ad hoc kinetic theory-type 
arguments. The present paper shows that similar predictions can be obtained by using a rigorous 
acoustic analogy approach together with actual measurements of the relevant acoustic source 
correlations. This puts the result on a firmer basis and enables its extension to new situations and 
to the prediction of sound at other observation angles. 

Introduction 

The prediction of aircraft exhaust noise continues to be a fruitful area of research. 
Computational demands still preclude the use of full scale DNS (or even LES 9 ) at the high 
Reynolds numbers and complex geometries of practical interest and the current emphasis 
remains focused on developing acoustic-analogy type approaches in which the Navier-Stokes 
equations are rearranged into a form that separates out the linear terms and associates them with 
propagation effects that can then be determined as part of the calculation. The non-linear terms 
are treated as “known” source functions to be determined by modeling and, in more recent 
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approaches, parameterized with the parameters being detennined from a steady RANS 
calculation. The “base” flow (about which the linearization is carried out) is usually assumed to 
be parallel and the resulting equation is usually referred to as a Lilley’s 1 equation. 

Most of the early approaches, including the original MGB“ approach, which neglect 
variations in retarded time in an appropriate moving frame coordinate system and assume the 
resulting moving frame correlation tensor to be separable into the product of (usually Gaussian) 
temporal and spatial components, significantly over-predict both the high and low frequency 
roll-off of the 90° acoustic spectrum. 4 Tam and Auriault 5 achieved much greater success in this 
endeavor by using ad hoc kinetic theory type arguments to derive their acoustic equation. 
However, Morris and Farassat 4 later showed that their result is equivalent to the usual acoustic 
analogy approach (at least in so far as its predictions of the 90° spectrum are concerned) with the 
primary difference being in the modeling of the source term, which amounts to assuming a (non- 
separable) functional form for the second convective derivative of the turbulence correlation that 
is similar to the functional form used for the turbulence correlation itself in the usual acoustic 
analogy approaches (but see refs. 30 and 31). Unfortunately, this forces the resulting turbulence 
spectrum to be singular and therefore incompatible with experimental observations, which would 
certainly be undesirable in any physics-based theory. 

The main purpose of this paper is to show that the Tam and Auriault 5 spectrum can be 
recovered from an acoustic analogy approach that uses experimentally based (i.e., physically 
realizable) source modeling together with an appropriate acoustic analogy equation. Reference 6 
shows that the Navier-Stokes equations can always be rearranged into the fonn of the linearized 
Navier-Stokes equations about a very general “base” flow, but with non-linear dependent 
variables, with the viscous stresses replaced by a generalized Reynolds stress and with the heat 
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flux vector replaced by a generalized stagnation enthalpy flux (which would be treated as 
“known” source strengths when the generalized equations are used as the basis for an acoustic 
analogy approach). This result leads to a fonn of Lilley’s equation with a modified source term 
when the “base” flow is taken to be a unidirectional transversely sheared mean flow, i.e., a 
parallel flow. 

It is generally agreed that the enthalpy fluxes, which correspond to the isentropic part of 
the pressure-density source in the Lighthill approach, are only important for hot jets ’ ’ ’ 
except, perhaps, at small angles to the downstream jet axis. 10 They are therefore neglected in the 
present analysis. The resulting equation is solved using a more or less conventional Greens’ 
function approach and the solution is then used to calculate the far field acoustic spectrum. 
Neglect of the enthalpy and viscous sources are the only approximations introduced at this stage 
of the analysis, but the results are then simplified by first neglecting variations in retarded time 
(in an appropriate moving coordinate system 11 ) and then introducing an axisymmetric turbulence 
model. “ ’ ’ Finally, an empirical turbulent source spectrum, based on recent measurements of 

17 . 

Harper-Bourne (in a low Mach number jet) is incorporated into the result and it is shown that 
the predicted far field acoustic spectrum at 90° to the jet axis is essentially the same as the one 
proposed by Tam and Auriault. 5 

The Acoustic Analogy Equation and its Far-Field Solution 
As indicated in the Introduction, reference 6 shows that the Navier-Stokes equations can 
be rewritten (for an ideal gas) as the linearized Navier-Stokes equations about a very general 
“base flow” but with different (in general non-linear) dependent variables, with the heat flux 
vector replaced by a generalized enthalpy flux and with the viscous stresses replaced by a 
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generalized Reynolds stress. This is a true acoustic analogy (in the Lighthill 18 ' 19 sense) in that it 
shows that there is an exact analogy between the flow fluctuations in any real flow and the linear 
fluctuations about a very general “ base flow ” due to an externally imposed “viscous ” stress and 
“heat flux ” vector. When the “base” flow is taken to be the unidirectional transversely sheared 
mean flow 

v. = 5 rl ty(x 2 ,x 3 ),p = p(x 2 ,x 3 ), 

p = p = constant (1) 


where x = {x, ,x 2 ,x 3 } is a Cartesian coordinate system, v = { v, , v 2 , v 3 } denotes the velocity,/? the 
pressure and p the density, the general equations reduce to the modified Lilley’s 1 equation 
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is the variable-density Pridmore-Brown operator 
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is the square of the mean- flow sound speed, y = the specific heat ratio, t denotes the time, 
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denotes the convective derivative based on 


(5) 


U, 


4 



is a generalized pressure fluctuation 
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is the generalized stress tensor and 
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is the generalized stagnation enthalpy flux. 
Here, 
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denote fluctuating quantities with h being the enthalpy and of and q' being the fluctuating 

viscous stress and heat flux vector, respectively, which are believed to play a negligible direct 
role in the sound generation process 18 ' 19 and are therefore neglected in the following. The 

fluctuating enthalpy flux r \' will also be neglected since, as noted in the Introduction, this 
quantity is generally considered to be unimportant for cold jets 1,7 ' 8 — except perhaps at small 
angles to the downstream axis. 10 

Then eq. (2) can be formally solved in terms of the free space Greens’ function 21 


G(jc,t| j,t), which satisfies 

LG(x,t\y,T) = b(x-y)8(t-T) (12) 



and has outgoing wave behavior at infinity, to obtain the following expression 

oo 

P 2 {x,t 0 )= J u (x\y;r\,t 0 +z 0 ) 

— OO V 

xR ijki ( j; h A )dydx\dx 0 (13) 

for the pressure autocovariance 22 (notice that p' reduces to // in the far field) 
p 2 {x,t 0 ) = —\ p' e (x,t)p' e (x,t+t 0 )dt ( 14 ) 

where V denotes integration over all space, T denotes some large but finite time interval, the 
propagation factor (x,*| y,x) is defined in appendix A and 

1 T 

R m{y^o) = ^ { P vyj{y,T) 

pv' k v'(y + x\,T + T o )dT ( 15 ) 

is the density-weighted, fourth-order, two-point, time-delayed fluctuating velocity correlation- 
with the indicated arguments referring to all three terms preceding the parentheses. The details 
are given in appendix A. 

Lighthill 1 8 ’ 1 9 pointed out that acoustic predictions tend to be less dependent on the 
detailed structure of the turbulence when variations in retarded time across the correlation 
volume A i) can be neglected, which is a reasonable approximation in a reference frame 

% = T[-\U C T 0 ( 16 ) 

moving with the local convection velocity, U c (y) ,of the turbulence. Ffowcs Williams 11 showed 
that this idea is best implemented by introducing the moving frame correlation tensor 

Rm ( r, 5. ■ ) = R m (r,i + > u c z 0 . r„ ) (r 7) 
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into the relevant pressure autocovariance fonnula (eq. (13) in our case) to obtain 


P 2 {x,t a )= J §i m {x\y\\ + \U c TQ,t 0 +T 0 ) 

— oo 

R tfki(rX*o) d y d l d * 0 ( 18 ) 

Our interest is in the far field spectrum 

4 ( • *) = J e m ° P 2 {x,t 0 ) dt 0 ( 1 9 ) 

271 J 
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which can, in principle, be calculated by taking the Fourier transform of eq. (18) and using the 
convolution theorem. 21 Unfortunately, this cannot be done in practice because the Fourier 
transform of R jjkl ( y;i /, z 0 ) , and therefore of R^ kl (y;£,, z 0 ) , does not exist 16 (see p. 179). This 

difficulty can be overcome by replacing R jjkl ( y;tj , z Q ) with 

Ryu (j;0,0)f? w (j + i/;0,0) where R fj (j;//,x 0 ) is defined by eq. (30). This does 

not change the radiated sound because R fj (p;0,0) R kI ( y + 0, 0) is steady and, therefore 

produces no sound. It then follows that 
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is the Fourier transform of y ;/ (we use capital letters to denote Fourier transform of the 
corresponding lower case quantity) and we introduced I a (x | y) ,the acoustic spectrum at x due 
to a unit volume of turbulence aty, i.e., 

7 «(*) = (22) 

v 

22 

in order to simplify the formulas. The relevant far field expansion of T ( y is given in 
appendix B. 

The only approximation made up to this point is the neglect of the enthalpy and viscous 
source tenns, but eq. (20) will depend on the turbulent source correlations only through 
= (23) 

V 

if variations in retarded time across the correlation volume are neglected, i.e., if 

i^ x y +g+iU c x; coj is assumed to be constant over the correlation volume. 11 However, the 

definition (17) implies that the integration variable in eq. (23) can be changed back to r| , which 
means that 

%jU ( y^o ) = I R rn {y>rj,%)dti (24) 

v 

i.e., the source correlation can be expressed in either the fixed or moving frame once the retarded 
time is neglected. (But introduction of the moving frame is necessary in order to get the correct 

t 0 -dependence in T * kl .) 

Equation (20) can now be written more simply as 

~ sin0F ^(*W F ^ 

(y;co(l-M c cos9))asx->°o (25) 


8 



where 


< 26 ) 

— oo 

is the spectral tensor of the source correlation and 



is the convective Mach number of the turbulence. This result shows that it is only necessary to 
model the overall spectral tensor itself and not the detailed two-point time delayed correlations of 
the turbulence. However, the radiated sound should still be relatively insensitive to the detailed 
turbulence structure even when the latter quantities are modeled (as is at least partially done 

below). This would not be the case if the moving frame had not been introduced before 

1 1 22 

neglecting the retarded time variations. ’ 

Our interest here is in the spectrum at 90° to the jet axis where cos0 = 0. Appendix B shows 

that 
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when Cq =c^ = constant , i.e., in the isothennal case. 
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The Quasi-Normal and Axisymmetric Turbulence Approximations 
To proceed further, we need to know something about the source spectral tensor The usual 
approach ’ ’ is to begin by assuming that the turbulence is quasi-normal (see ref. 10) in order 
to obtain some relations among its components. It then follows that (see comments preceding 
eq. (20)). 

R m ( v; n, x 0 ) = R ik (y; % x 0 ) R jt (y; % x 0 ) 

+ R n ( T; % \ ) R jk (y;r\, x 0 ) (29) 

where 

1 T 

R ij {yw 0 ) = ^ J Vp v 't (ta)Vpq 

{y + %x + T 0 )di (30) 

is the second order correlation. To further reduce the number of independent components it is 
usual to assume some kinematically possible symmetric form for the second order correlations. 
Early studies assumed the turbulence to be isotropic, but that turns out to be incompatible with 
the Harper-Bourne 17 measurements that will be introduced below. The simplest assumption 
compatible with his results is the one introduced in references 12 and 13, namely that the 
turbulence is axisymmetric which implies that 16 

) = MKj + B A) + c„s„s, j 

+A(M; + M.) (31) 

where An, Bo, Co, and D 0 are functions ofy, r 0 , and r\ ± ; Ao, B 0 and Co are even functions r\ ± and 

Do is an odd function of this quantity. This model is chosen because it is the most general of 
those whose mathematical properties have been studied in the literature and because it reflects 
the fact that the cross flow velocity components tend to be much more similar to one another 
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e ° { (^12 + ^1 1^22 )^n dx 0 
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( 33 d) 


are seemingly independent spectral functions. However, the coefficients Ao, Bo, Co, and D 0 are 
not all independent and when compressibility effects are neglected (i.e., when p is treated as a 
constant) these turbulence correlations can be expressed in terms of two independent scalar 
functions of y, x 0 , rp , and p, , say a and b, both of which are even functions of the latter 
variable . 14 ' 13 ’ 24 The resulting expressions for the two point correlations are given in appendix C, 
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where they are used to express the integrals over V in eq. (33) in terms of a and b. Equations 
(C — 1) to (C-4) suggest that these latter quantities will scale like 

b = pu^B(f\ 1 ,f\ l )l? L /2 (34) 

and 

s ~ a ~\y\ l = P u 2 D ( f \±’ f \i) ( 35 ) 

where 


ill = VA 

(36) 

ill =il±/A . 

(37) 


I, and L ± denote characteristic stream-wise and transverse length scales of the turbulence, B and 
D are 0(1) functions of the indicated arguments, 

pif=R n {0, 0) (38) 

and 

pu\ = R 22 (0,0) (39) 

Turbulence measurements suggest that 

8 - —Jfe (40) 

4 A 

17 

ought to be small. In fact, Harper-Bourne’s measurements (to be discussed below) suggest that 
e- 2.7xl0“ 2 . The scalings (34) through (40) are inserted into eqs. (C-5) through (C-9) of 
appendix C where it is shown that 
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Lacking any specific data to the contrary, it seems reasonable to suppose that 
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The Harper-Bourne Spectrum 

The results can not be made more explicit without inputting more specific information 
about the turbulence structure. This is accomplished with the aid of some recent measurements 17 
of the two point fourth order stream-wise velocity correlation spectra along the centerline of the 
mixing layer in a low Mach number jet, which would most closely correspond to 

H o {y,ri,<a) S -]e- ,m -R^(y,%x o )dx 0 (47) 

K J 

— OO 

with the quasi-normal approximation that is being used in the present analysis. 

Harper-Bourne 17 divided H 0 into the three components (see his eq. (2.5) and (2.7) on 

p. 2) 
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H 0 = H 0 (j,0,co)i? 
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where l x , l L are the spectral stream- wise and transverse length scales (not necessarily the same 


as the time domain length scales L { and L ± introduced above) and 
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p 



(49) 


No assumption is made about the decomposition of the correlations into products of their space 
and time components with this approach. 

The first factor can be evaluated from his measurements of f? llu (y,0,x 0 ) , which are well 


represented by the exponential 



-implying that 


Ho{y,0,o>) 
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Inserting these into eq. (43) and using the result in eq. (44) shows that 
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and /, and /_ are defined by 
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( 53 ) 
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Harper-Bourne’s measurements seem to indicate that j,r|p(o) behaves like e , which, as 

pointed out by Khavaran et ah, is inconsistent with Hinze’s” parabolic vertex requirement. We 
therefore apply a Gaussian filter to his measured value of R, say R , to obtain 






Vne^ 1 ' erfc.f ^ 


which smoothes out the cusp at ru =0 but, otherwise behaves similarly to R m (note the data 


scatter around = 0 in Harper-Bourne’s fig. 14(a)). The Error function factor, e • erfc.f 


has been inserted in the denominator in order to insure that R still satisfies the normalization 
condition, f?(j,0,0) = 1, when R m fj,rjpCA = e ~^ . In fact it is easy to show more generally 


that, with this choice of R m f j,T|p0j , 
e~^' erfc.f ^ ^ 1 +e^ erfc.f ^ ^ 
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which behaves like e n ' for large r\ { but has continuous slope at r\ { = 0 . 

Inserting eq. (55) into eq. (52), using the convolution theorem, and inserting the result 
into eq. (51) yields 
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which reduces to eq. (51) when (3 = 0. 


Comparison with the Tam and Auriault result 

The fonn 

R m = e -^e^ % 1 (59) 

agrees fairly well with Harper-Bourne’s measurements when = 0 and T|_ 0 and also 

whenr| ± = 0 and Hj ^ 0. Inserting this into eq. (58), integrating over T| and inserting the result 
into eq. (57) shows that 
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The form (59) is consistent with the assumption that f? / 2 = D = a (r[ 1 , x o ) e~^ 2 when /, and / ± 
are taken to be constants (i.e., independent of co), which then implies that the ratio of integrals in 
eq. (46) is equal to Vi and that T 2 = 3 in eq. (43). So that for r = Vi, y = 1.4, and U c = 0.65 Uj 


M 1 

2 Uj 2 l s Dj 


(61) 


where / s ~ 0.345 denotes the characteristic length scale of the mean velocity divided by the jet 
diameter D , . Harper-Bourne’s measurements show that / ± and 1 1 are relatively constant and that 
the latter can be approximated by /, = .69 ID ; at sufficiently low frequencies. It therefore follows 


that the square bracket in eq. (60) is approximately equal to unity, which means that the resulting 
90° acoustic spectrum should be of the fonn 



e-l^^co 2 

(to 2 +?t 2 ) 


(62) 
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which is the same as the spectral form proposed by Tam and Auriault, who take — - = —^^~ 

U 2e u 


and X = y c ^ (see eq. (53) of ref. 4) where c l and c\ adjustable constants and k and £ (not to be 
confused with the £ in eq. (40)) are determined from a k-z RANS calculation. Notice that this 

result behaves like or as co— > 0 and converges like e _, P 0)/ N) as ^ ^ 

However, Harper-Bourne demonstrates that the separable form (59) does not work at 
oblique separations where and T|_ are both non-zero and, more importantly, his fig. 13 shows 

that /j and l L are only constant at relatively low frequencies with the scaled length scales ( and 


l L becoming constant as cg-^ °° . It, therefore, follows from eq. (57) and the asymptotic behavior 


of the Error function (ref. 29, p. 298 #7.1.23) that / 0) (*|j) converges like 


e -(N/r #) 2 


when 


K ~ const, and co — > °o s i.e., it has the same high frequency behavior as Tam and Auriault result 
(62). The general fonnula (57) therefore coincides with that result at both high and low to 
moderately-low frequencies-even with frequency dependent length scales. 

The two results differ at intermediate to moderately high frequencies by the factor 
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(where /[°^ = .69 ID 7 and /| 0 ^ are the initial constant values of /, and l ± ), which is no longer 
equal to unity at these frequencies but becomes constant as co — > °o . On the other hand, 
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Khavaran et al. 27 point out that the Tam and Auriault 5 data comparisons do not account for 
atmospheric attenuation. Figure 1 ( taken, in part, from ref. 27) shows that the latter comparisons 
(curve 2) under predict the high frequency portion of the spectrum for a Mach number 0.5 cold 
jet when atmospheric absorption effects are included. But multiplying the result shown in 
curve 2 (which involves a summation over the jet) by the overall (or average) correction factor 
(63) (with /, and/_ given by the fonnulas in Harper-Bourne’s fig. 13), which is hopefully 

relatively independent of the choice of source point, leads to curve 1, which is in much better 
agreement with the data. It would, of course, be better to sum the actual formula (57) over the jet 
with appropriate local values for the parameters, but since the measurements are all at a single 
location, this does not seem to be warranted at the present time. 

Harper-Bourne obtains the best fit to his data with the non-separable form 

R m = e~^ (64) 

which can be inserted into eq. (57) via eq. (58) to obtain 

Xp (fe ( ' |i “ ‘ * L- + (u k)-1 

C y 

j j (V + 0) )co[l + (2jtTf J 

xtu/77 2 (65) 

which has the same high and low frequency behavior as eq. (60) but is slightly different from it 
at intermediate frequencies. More in depth data comparisons may be needed to determine if this 
result will provide a better representation of the data. 
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Discussion 


As noted in the Introduction, Morris and Farassat 4 argue that, for the same mean flow 
calculation, the Tam and Auriault result (62) yields a better prediction of the 90° acoustic 
spectrum than methods based on the usual acoustic analogy approximations. They also show that 
this type of spectrum can be obtained within the acoustic analogy frame work if the usual source 

D 2 

modeling is applied to — - R jjkl rather than to R jjkl as was done above. But this would imply 

Dx~ 


that the R ijkl spectrum becomes infinite like co " as ©->0, which is certainly not consistent with 


experimental results. The conventional acoustic analogy approach leads to a multiplicative factor 
of co 4 in the spectrum 22 (rather than the to 2 factor in eq. (62)) due to the quadrupole nature of the 
source). The co 2 factor is characteristic of a dipole source. In the present approach, the energy 


dU 


in- 


equation introduces the dipole type source tenn (y - 1) — -e Xi in the acoustic analogy 

dx, Dr 


eq. (2) which combines with the quadrupole source and a portion of the physically realizable 
spectral function (eq. (48)) that arises from the non-separable nature of the source to produce the 
correct spectral form. 

It is important to keep in mind that while the 90° spectrum is useful from a diagnostics 
point of view, the maximum acoustic intensity occurs at relatively small angles to the 
downstream axis for most (if not all) high speed air jets. The predictions of the present model 
would differ from the Tam and Auriault 5 results at these angles or, for that matter, from any 
modification of the conventional acoustic analogy. 
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Appendix A 

The Pressure Autocovariance 

The formal Green’s function solution to eq. (2) can be written as 


oo 

Pe (*>0 = J J P v i v 'j Y,y (x,t\y,r)dydr 


V — OO 

where pv'v' is evaluated at y,r , 
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\^y m a T« Dx d y„ 

du D 2 


dv, Dx 2 


G(xy|j,x) 


dU d 
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where all the terms in square brackets are part of the differential operator and we have integrated 
by parts to transfer the derivatives from the source term to G. Inserting this into eq. (14), 
changing integration variables to t\ = t - x 2 and 
Xj = x 2 -x , and introducing eq. (15) yields 


1 i ~ 

P 2 = — J \\\\ln( x W + t o - x )YA/(^|T^- x 2 )p v Ay (T^)pv>;(ji,x 2 )JjJji dxdx 2 dt 

^ -T— V 

-M Y ij +t Q +Xj) y kl (x|jj ,tx)R ijkl - y,^)dy dy x dt x dx x (A -3) 
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which upon introducing the separation vector 

*1 = yi~y (A -4) 

can be written more compactly as eq. (13) with 
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yijki{x\y;ijA +^o) 

oo 

= J y i j{x,y,ti+t 0 +^ 0 )y kl (x\y+n,t l )dt l 

— oo 

(A -5) 

The quantities y l , t x , x, are dummy integration variables and we are using the fact that G and 
therefore y (/ depend on t and x only in the combination t - x . 
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Appendix B 

Far Field Expansion and Neglect of Retarded Time 
Taking Fourier transforms of eq. (12) with respect to (x, -y x ) and (t-x) shows (since G 
can depend on these quantities only in these combinations) that 

H x ±-y±) 


A Go 
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where 
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is the reduced Rayleigh operator and 
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(B-3) 

where x ± and y ± have the obvious meaning. Inserting this into eqs. (21) and (A-2) shows (after 
some rearrangement) that 
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? ? Y 1 ? 

^in ^ jm 2 y nm 
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(< o - a - i /(^)) 2 9y. 


2 ?-e'^G 0 M-il> u (y-l) 


du 7 e ik{x '- y ' ] G. 
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L(®-kU(y L )) 


°-xdk (B — 4) 
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This result is much simpler in the far field where U (x ± )—>(), c/ — > c 2 = constant, and 




e -x±J/c 2 -(w/c „) 2 





as jc_l — > °° (B - 5) 


Here, x ± = IxJ and (p = tan 1 — denotes the circumferential angle. Inserting this into eq. (B-4) 

.v, 


and using stationary phase to evaluate the integrals now shows that 

cose)/c„ l 27l/0)sin9 


T ij 
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iaix/c^ 


2ra'cosin0 


r A x M = 


( y — 1 ^ 

2 2 / 1 s ? 

u jm ^ u ij ^ nm 

\ Z J 


X- 


c 2 /of 


dy m [m(j_l)cos0-1] 

/5i ; (y-i) df/ i 

co 3v y [M(j ± )cos0-l] 


Xe -m,v, cos9/c„ ~ 


cp, — cos©,®!^ 


as x — > °° (B - 6) 


where x = Lvl is the radial coordinate, 0 = arcsin 




is the polar angle of the observation 


point and 

M = U( y ± )/ Coo (B — 7) 

is the “acoustic” Mach number at the source point. 

Since eq. (B— 1) reduces to Helmholtz’s equation when k = 0, G a (jc ± |j ± ;0,co) is given 
by 21 (seep. 811) 
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when cos0 = 0 and c 2 0 =c^ = constant , it follows from eqs. (B-5), (B-6), and (25) that 
do{(p,0,a\y ± ) 



and therefore that eq. (28) holds. 
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Appendix C 


Reduction of Correlation Coefficients 


It is shown in references 15, 24, and 14 that 


R 2 9 — a H — — 


r\l da d 2 b 

hi dr\ ± dri. 
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r \ 2 ^ 2 b 

hi ^hi^hi 


(C-3) 


D 1 d db 
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so that integration with respect to the circumferential coordinate shows that 
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Inserting eqs. (34) and (35) into these results and integrating over ri_|_ shows that 


Ji? 2 2 clx\ 3 j /? 2 2 3 dx\ 


2k L x I} l ( p wj 2 ) 2^11 ( p» 2 


+ 0(e J ) 


7 7f„ 3D ' 

i!KJ 

l 

f\ L df\ L d% + o(e 2 ) 


(C-10) 

J R n dr\ 


V 

„ / — 1\ 2 

=0(e 2 ) (C — 11) 


2kL { L\ pwf 


J dr\ 


2 2 


oo oo 

= -{ J + o(e 2 ) 

"^»o 


(C-12) 


2kL x I} l pw 2 


=§//" 


3 ^ 


— OO 0 


3f) 

where we have put 


+6>(e 2 ) (C-13) 


B = - 
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ill df\ A 
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drij 
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and r is defined by eq. (42). 
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